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Abstract 

The X-ray emission from hot gas in bright elliptical 
galaxies often extends far beyond the radius of the 
stellar system. This "circumgalactic" gas accounts for 
most or all of the large spread in X-ray luminosity 
among ellipticals having similar optical luminosities 
Lb. 

We have developed gas dynamical models describ- 
ing the evolution of gas within and around ellipticals 
beginning with an overdensity perturbation in a sim- 
ple flat cosmology. At some early time we form the 
stellar galaxy and release supernova energy, conserv- 
ing dark and baryonic matter. We follow the subse- 
quent evolution of intergalactic and interstellar gas to 
the present time. These models confirm that hot gas 
density and temperature distributions currently ob- 
served in massive, group-dominant ellipticals can be 
understood as a combination of intergalactic gas that 
has flowed into the galaxy group over time with gas 
lost from galactic stars. Furthermore, if the hot gas 
and dark matter halos are subject to differential tidal 
truncations or mass exchanges between group mem- 
bers, then the observed correlation between L^/Lb 
and the relative sizes of galactic X-ray images can be 
understood. The distribution and physical properties 
of hot interstellar gas observed in massive ellipticals 
today are sensitive to the cosmic baryon fraction, the 
time of maximum star formation and the amount of 
"feedback" energy delivered to the gas by Type II 
supernovae at the epoch of galaxy formation. 
Subject headings: galaxies: elliptical and lenticular - 
galaxies: formation - galaxies: evolution - galaxies: cool- 
ing flows - X-rays: galaxies 

For many years X-ray observers have noticed that 



the X-ray luminosities of elliptical galaxies span a 
huge range for galaxies of similar optical luminos- 
ity Lb (Fabbiano 1989; Eskridge, Fabbiano, & Kim 
1995). Attempts to find correlations between L^/Lb 
and other galactic parameters had been largely unsuc- 
cessful. However, in a recent compilation of ROSAT 
data, Mathews & Brighenti (1998) discovered that 
Lx/Lb correlates strongly with the size of the X-ray 
image v^x /^e i where t^x 

is the projected radius that 
contains half the total Lx and is the optical half- 
light radius. A correlation between Lx/Lb and r^x/re 
was suspected because the hot gas in well-resolved, X- 
ray luminous ellipticals often extends far beyond the 
optical images where it could be rather stochastically 
influenced by tidal mass exchanges. 

The very extended X-ray halos around many bright 
ellipticals also indicate that most of this gas is un- 
likely to have been ejected by old galactic stars during 
their post-main sequence evolution. Furthermore, the 
average gas temperature in bright ellipticals exceeds 
the equivalent stellar temperature by 1.5 (Davis 
& White 1996), suggesting again that the hot gas is 
in virial equilibrium in the more massive and more 
extended dark halos surrounding the galaxies. Using 
detailed hydrodynamic models, Brighenti & Mathews 
(1998) showed that the extended hot gas density and 
temperature profiles in X-ray bright ellipticals can be 
understood only if most of the hot interstellar gas is 
very old, dating from the epoch of galaxy formation, 
and has not come solely from mass loss from galactic 
stars. 

In this letter we combine these two lines of inquiry. 
First, using a simple cosmological model we show that 
the additional "circumgalactic" gas required to ac- 
count for large values of Lx/Lb can result quite nat- 
urally from secondary infalling intergalactic gas that 
has accumulated in the outer dark halo over a Hub- 
ble time. Second, we demonstrate that the correlation 
between Lx/Lb and r^x/re can indeed be generated 
by tidal truncations of these same models. 

For a fully self-consistent model of the global evo- 
lution of gas in ellipticals, it is necessary to begin 
with a galaxy-group sized perturbation in an expand- 
ing universe, allow for the formation of the stellar 
elliptical, and include the release of supernova energy 
that accompanies early star formation. For simplicity 
we consider here a localized density perturbation in a 
simple flat cosmology with various baryonic fractions 
flf)/Q. After a few gigayears, when enough baryonic 
matter has accumulated to create the giant ellipti- 
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cal, wc construct the stellar galaxy with a dc Vau- 
couleurs profile and simultaneously release the energy 
from Type II supernovae. By forming the galaxy in 
this instantaneous fashion, we circumvent the dynam- 
ically complex merging processes that must have oc- 
curred (e. g. Merritt 1985). After the stellar system 
is formed, baryonic and dark matter continue to flow 
into the galactic perturbation. The baryonic mat- 
ter shocks and compresses to approximately the virial 
temperature in the halo, T ^ 10^ K, similar to inter- 
stellar gas temperatures observed today. With this 
simple model we conserve baryonic and dark matter 
and allow properly for the influence of supernova en- 
ergy. We study the subsequent gas dynamical evolu- 
tion of the intergalactic and interstellar gas in detail, 
allowing for stellar mass loss and Type la supernovae. 
The objective in this calculation is to reproduce the 
density and temperature profiles observed in the hot 
gas within and around bright ellipticals today. A 
more detailed discussion of our hydrodynamical mod- 
els has been submitted to the Astrophysical Journal. 

At the present time only about a dozen bright, 
nearby ellipticals have been spatially resolved in soft 
X-rays. For these galaxies, the X-ray surface bright- 
ness and projected spectral variation can be inverted 
to determine the electron density n(r) and temper- 
ature T(r) profiles with physical radius. Assum- 
ing approximate hydrostatic equilibrium, the density 
and temperature gradients can be used to determine 
the total internal mass M(r) (Brighenti & Mathews 
1997). As a guide in evaluating the success of the 
models we discuss here, we choose NGC 4472, a typ- 
ical bright E2 elliptical in Virgo that is the dominant 
elliptical in a local galaxy subgroup. Data for NGC 
4472 is available from Einstein HRI (Trinchieri, Fab- 
biano, & Canizares 1986) and ROSATSRl and PSPC 
(Irwin & Sarazin 1996). The interstellar gas density 
profile for NGC 4472 has been determined out to al- 
most 17re « 150 kpc (at a distance of 17 Mpc). NGC 
4472 appears to be colliding with a small dwarf galaxy 
and interacting with more extended Virgo intercluster 
gas (Irwin & Sarazin 1996; Irwin, Frayer, & Sarazin 
1997). This may explain the azimuthal asymmetry in 
the X-ray image beyond ^ 3.5 kpc. In spite of this 
undesirable asymmetry, in these times of pre-AXAF 
resolution NGC 4472 is one of the few normal bright 
ellipticals in which the stellar component is clearly 
seen in the mass profile determined from X-ray ob- 
servations; this is useful for our models to accurately 



evaluate the stellar potential, the gas contributed by 
stellar mass loss, etc. The X-ray profile in NGC 4649 
is also sensitive to the stellar potential and its X- 
ray properties are very similar to those of NGC 4472 
(Brighenti & Mathews 1997), but X-ray observations 
in NGC 4649 are not as spatially extensive as in NGC 
4472. On balance, we adopt the global, azimuthally 
averaged X-ray properties of NGC 4472 as generally 
representative of other X-ray bright ellipticals. 

The total mass distribution for NGC 4472 found 
by Brighenti & Mathews (1997) can be fit quite well 
with a de Vaucouleurs stellar core out to ^ and an 
Navarro-Frenk- White halo (NFW) beyond (Navarro, 
Frenk, & White 1996). In making these fits we use 
a total stellar mass of = 7.26 x 10^^ Mg which 
is based on Lb = 7.89 x 10"'^'' ^b,0 f^nd the mass to 
light ratio M^/Lb = 9.20 of van der Marel (1991); 
this M^/Lb is in fact verified by the X-ray image of 
NGC 4472 (Brighenti & Mathews 1997). For the one- 
parameter NFW halo we use a virial mass of Mh = 4 x 
10^^ Mq which corresponds to an NFW concentration 
parameter of c = 10.47. The total mass distribution, 
stars plus dark halo, agrees very well with the X-ray 
determined M(r), but is slightly too massive near ~ 
Te suggesting that the dark halo may be less centrally 
peaked than NFW predict. 

In our hydrodynamical model for the evolution 
of the intergalactic and interstellar gas, we adopt a 
very simple flat cosmology with = 1 and 7? = 50 
km s~^ Mpc~^, corresponding to a current universal 
age of i„ = 13 Gyrs. Primordial nucleosynthesis re- 
stricts the baryonic mass fraction to a narrow range, 
nb{H/bQf « 0.05±0.01 (Walker et al. 1991). We as- 
sume that both baryonic and dark matter converge to- 
ward the proto-4472 galaxy perturbation as described 
by the self-similar solution of Bertschinger (1985). In 
this model all flow parameters depend on a single sim- 
ilarity variable A = r/rta{t) where rta{t) oc t^/^ is the 
turnaround radius where the cosmic flow velocity van- 
ishes at time t. 

Before we form the galaxy at time the cold bary- 
onic gas within rta{t) flows exactly with the dark mat- 
ter until it encounters an accretion shock. The time- 
dependent dark matter potential is described by an 
outer Bertschinger flow, similar to that into a cen- 
tral singularity, and an inner stationary NFW halo 
characterized by parameters that fit M (r) for NGC 
4472. At each time the NFW and Bertschinger mass 
distributions are matched at some radius rc{t) that 
conserves the total mass of dark matter. Outside 
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both Tc{t) and the ac;crc;tion shock, the mass flux 
in the baryonic component at fixed radius varies as 
pu oc t~^/^ so most of the halo gas is accumulated 
at early times. The accretion shock forms at a very 
large radius, Vgh ^ 1 Mpc for t > 6 Gyrs. Within the 
accretion shock a cooling flow is is established and the 
gas velocity is very subsonic. 

At some time when a sufiicient number of baryons 
has collected in the center of the flow, we form the de 
Vaucouleurs stellar mass distribution (r) matching 
NGC 4472. This is done by reducing the gas density 
within the accretion shock radius rsh{t*), lowering its 
mass by M*. Also at this time the remaining gas in 
r < rsh{t*) is heated by the Type II supernova en- 
ergy from massive stars, m > 8 Mq. We assume a 
Salpeter power law IMF from mi = 0.08 to to„ = 100 
Mq for which rju = 6.81 x 10~^ SNII are produced per 
Mq of stars formed, each of energy Egn = 10^^ ergs. 
The amount of SNII energy received by the hot gas 
depends both on the uncertain IMF-dependent value 
r]ii and the efficiency that this energy is transferred 
to the hot gas and not lost by radiation. For this rea- 
son we adopt ijstdEsnM^ with rjstd = 6.81 x 10~^ as 
a reference energy which can be adjusted as needed. 

After time the stars lose mass as described by 
a^{t) = 4.7 X 10-20[t/(t^)]-s where i„ = 13 Gyrs 
is the current time, tg = tn — and s = 1.26. We 
assume that the gas is heated after by Type la 
supernova, each of energy Esn, at a rate SNu(i) = 
SNu{tn){t/tn)~P with p = 1 and SNu(i„) = 0.03 su- 
pernovae per 100 years per 10^*^X5,0. 

The gas dynamical equations we solve are 
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where s = 3kT/2pmp is the specific thermal energy; 
p, = 0.62 is the molecular weight and nip the proton 
mass. Mtot{r,t) is the total mass of stars, dark mat- 
ter and hot gas within radius r. The mean injection 
energy of gas from stars is = 3kTo/2pmp where 
To = (a*r* + asnTsn)/a and a = a* + A(r) 
is the coefficient for optically thin radiative cooling. 
These equations are described more fully by Brighenti 



& Mathews (1998). We solve them numerically using 
a spherical ID Eulerian hydrocode (ZEUS) with log- 
arithmic spatial zoning. 

In Figure 1 we illustrate the results of several gas- 
dynamical calculations at time tn = 13 Gyrs using 

= 2 Gyrs, = r^std and fib = 0.05 or 0.06. 
Considering the relative simplicity of our model, the 
agreement with both the temperature and density dis- 
tributions observed in NGC 4472 is excellent. The 
positive gas temperature gradient within r ^ 3 — 5re 
kpc observed in NGC 4472 and other bright ellipti- 
cals (Brighenti & Mathews 1997) occurs as hot in- 
flowing halo gas mixes with the slightly cooler gas 
at ~ T* that characterizes gas ejected from stars or- 
biting within a few r^- The apparent sensitivity of 
our results to fif, is of considerable interest. However, 
we emphasize that the influence of the parameters 
fife, T]!! and is degenerate. For example, solutions 
with decreased by ^ 0.01 are similar to those with 
Tjii increased by ~ 2 or with decreased by ~ 1 
Gyr. This degeneracy can be removed with more 
detailed models that include the abundance of iron 
and other elements produced in supernovae (Renzini 
1998). Nevertheless, the Ob = solution shown in 
Figure 1, for which the only source of gas since is 
stellar mass loss, is clearly inadequate. We conclude 
that cosmic inflow of gas is essential in understanding 
the properties of hot gas observed in giant ellipticals 
today. Moreover, this extended "circumgalactic" gas 
is sensitive to important cosmological parameters and 
to the details of SNII energy release during the epoch 
of galaxy formation; solutions with little or no SNII 
energy provide very poor fits in Figure 1. 

In Figure 2 we show the current distribution of the 
baryon fraction = Qb/fi with radius for the fib = 
0.05, T]!! = r]std solution. The baryonic cavitation 
around the galaxy is a relic of the SNII energy released 
at when gas was pushed out of the central parts 
of the flow. This suggests that determinations of fb 
based on X-ray observations in the range ^ 100 — 300 
kpc could either underestimate or overestimate the 
true cosmic value of fib- Obviously, the variation of 
/b(r) in r > 300 kpc shown in Figure 2 cannot be 
relevant to NGC 4472 because of the gravitational 
influence of the Virgo cluster. 

As a further confirmation of the overall success 
of our models, we compare in Figure 3 the locus 
of truncated models with observed ellipticals in the 
{Lx/Lb, rex/''e)-plane. The large circle, cross and 
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square at the right arc the current loci of untruncatcd 
models with = 0.04, 0.05 and 0.06; we have made 
these models somewhat more gas rich than NGC 4472 
by reducing ry//. When these models are truncated at 
time ttr = 9 Gyrs the current position of the galaxy at 
time tn = 13 Gyrs moves progressively along the ob- 
served sequence of ellipticals as the truncation radius 
rtr decreases from 500 kpc to 100 kpc. At the time 
of truncation, all gas and dark matter is removed at 
r > rtr, causing a rarefaction wave to propagate into 
the galactic ISM; the cosmic inflow is stopped after ttr 
and a new equilibrium is established in a sound cross- 
ing time ^ 1 Gyr. These results arc insensitive to 
the time ttr when the truncation occurred as long as 
ttr ^ 6 Gyrs. This insensitivity to ttr arises from the 
rather small changes in the gas flow variables within 
the accretion shock until the present time, the large 
radius of the accretion shock at early times {rgh ~ 400 
kpc at i = 4 Gyrs), and the declining rate of intcr- 
galactic gas accumulation with time. The density and 
temperature profiles at i„ for the truncation solution 
that best matches NGC 4472 (fib = 0.05, rt,. = 400 
kpc) are shown in Figure 1. This suggests that NGC 
4472 may have been tidally truncated in the past. 
Finally, we expect that purely gaseous truncations 
produced by ram pressure stripping in dense clusters 
would also reduce both L^/Lb and Vex/fe (Sakelliou 
& Merrifield 1998). 

Our main conclusions: 

• The extended "circumgalactic" hot gas surround- 
ing massive ellipticals can be understood as a natu- 
ral consequence of inflowing intergalactic gas over the 
Hubble time. This gas cannot have been produced by 
mass loss from stars within the elliptical itself. 

• Gas provided by stellar mass loss virializes to a 
relatively lower temperature ~ T^, and, when mixed 
with hotter gas flowing in from the outer halo, creates 
the positive gas temperature gradients commonly ob- 
served within ~ 3re (Brighenti & Mathews 1997). 

• Currently observed hot gas density and temperature 
profiles are sensitive to ^b, to the time of galaxy for- 
mation, and to the total "feedback" energy released 
by Type II supernovae as galactic stars formed. 

• Because of the SNII-driven outflow just after the 
time of galaxy formation, the cosmic baryonic frac- 
tion fife/il cannot be accurately measured from X-ray 
observations or dynamical arguments within several 
hundred kpc. 



• The fraction of hot gas currently within ^ r,. that 
comes from stellar mass loss is strongly dependent on 
Sib and the SNII energy. This must be considered in 
interpreting the abundance of iron and other elements 
in the hot gas. 

• The enormous range in L^/Lb among elliptical 
galaxies is due; in large part to differential truncation 
of the outer halos. The truncation could be caused by 
massive nearby group member galaxies or by another 
nearby group or cluster. It is likely that NGC 4472 
has been truncated, perhaps by nearby group com- 
panions or by the Virgo cluster. If the wide range in 
Lx/Lb were due to Type la supernova-driven galac- 
tic winds, as is often suggested, the L^/Lb, rgx/re 
correlation would be difficult to explain and the hot 
gas iron abundance would depend on L^/ Lb- 
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Fig. 1. — Gas density {Upper panel) and tempera- 
ture [Lower panel) of several models at i„ = 13 Gyrs 
compared to observations of NGC 4472. Solid line: 
fife = 0.05, = 2 Gyrs, r/// = r]std\ Dashed line: 
Vtb = 0.06, t* = 2 Gyrs, rjn = rjsid', Short dashed 
line: flj, = 0.05, = 2 Gyrs, rju = O.'Srjstd^ Ur = 9 
Gyrs, rtr = 400 kpc; Dash-dotted line: Cli, = 0.00, 
= 2 Gyrs, rju = 0. Filled circles are Einstein 
HRI gas density observations from Trinchieri, Fab- 
biano & Canizares (1986); open circles are ROSAT 
density and temperature observations from Irwin & 
Sarazin (1996). 

Fig. 2. — Solid line: Variation of the baryon fraction 
/b at time i„ = 13 Gyrs. Dashed line: The baryonic 
contribution of the gas alone. The model parameters 
are fib = 0.05, = 2 Gyrs, and ?7j/ = -qsui- At 
very large radii the baryonic fraction approaches fib = 
0.05. 

Fig. 3. — A plot of the (La; /Lb, rex/?'e)-planc show- 
ing the correlation among observed ellipticals taken 
from Mathews & Brighenti (1998). The large circle, 
cross and sqiiarc at the right arc loci of untruncatcd, 
gas-rich models at t„ with the following parameters: 
(Qfc, t*. Till) = (0.04, 3, 0.1), (0.05, 2, 0.3), and (0.06, 
2, 1) respectively. The smaller symbols that extend 
from the untruncated solution toward the lower left 
show the current {t = tn) loci of models truncated 
at time fj,. = 9 Gyrs at five decreasing tnmcation 
radii rtr = 500, 400, 300, 200, and 100 kpc. As rtr 
decreases the current locus in the [L^/Lb, Vex/'t'e)- 
plane moves from the untruncated solution almost 
exactly through the observed correlation. The filled 
square at the lower left is the locus of NGC 4472 if 
stellar mass loss creates all of the interstellar gas (the 
fife = solution). 
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